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Stochastic acceleration of charged particles due to their interactions with plasma waves 
may be responsible for producing superthermal particles in a variety of astrophysical 
systems. This process can be described as a diffusion process in the energy space with 
the Fokker-Planck equation. In this paper, a time-dependent numerical code is used 
to solve the reduced Fokker-Planck equation involving only time and energy variables 
with general forms of the diffusion coefficients. We also propose a self-similar model for 
particle acceleration in Sedov explosions and use the TeV SNR RX J1713. 7-3946 as an 
example to demonstrate the model characteristics. Markov Chain Monte Carlo method 
is utilized to constrain model parameters with observations. 
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1. Introduction 

Observations of high-energy emission from many astrophysical sources imply the 
presence of significant populations of relativistic particles. Relativistic cosmic rays 
(CRs) with energies lower than ~ 10 15 eV are commonly believed to be accelerated 
at the shock of supernova remnants (SNRs) in our Galaxy UJ Evidence of particle 
acceleration by shocks of SNRs first comes from radio observations of the remnants, 
where the radio emission is produced by relativistic electrons via the synchrotron 
process. The discovery of synchrotron X-ray emission from SNR 1006 reveals accel- 
eration of TeV electrons by these shocks^. However, direct evidence of proton and 
ion acceleration up to the spectral "knee" by these shocks remains elusive. 

Nonthermal distributions are naturally produced via the Fermi process. The 
first-order Fermi mechanism corresponds to a convection term in the energy space. 
The second-order process, as it was originally conceived by Fermi, involved the 
stochastic acceleration (SA) of particles scattering with randomly moving magne- 
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tized clouds Later refinements of this idea replaced the magnetized clouds with 
magnetohydrodynamic (MHD) waves SI. The second-order stochastic Fermi pro- 
cess recently finds applications in a wide range of astrophysical environments (See 
Refs. ISHII for more details). 

One commonly assumed agent of SA is plasma wave turbulence, which is ex- 
pected to be present in nonequilibrium conditions of highly magnetized plasmas. 
The wave-particle interactions in magnetized plasmas have influences on the dynam- 
ics of energetic charged particles and can play a significant role in the propagation 
and acceleration of these particles. Charged particles which spiral along magnetic 
field lines are accelerated through resonant interactions with plasma waves. This 
process can be described as a diffusion process in the energy space through the 
Fokker-Planck (F-P) equation with a diffusion coefficient whose magnitude and form 
depend on the power spectrum and other characteristics of the plasma turbulence. 

The resultant equation in general is complicated and one has to resort to some 
simplifying approximations. If we are mainly interested in the spatially integrated 
results, one may obtain a reduced F-P equation with time and energy as the vari- 
ables. There have been a few studies for analytic solutions of this reduced F-P 
equation in simple cases However, numerical methods need to be developed 
to solve more general forms of the reduced equation. Park & Petrosion (1996) 
studied in detail the numerical methods for the reduced F-P equation. In this paper, 
we use their f ully implicit Chang-Cooper method, which is proposed by Chang & 
Cooper (1970) OH. Our primary goal in this paper is to present a detailed deriva- 
tion of the exact time-dependent numerical solution for the re duce d F-P diffusion 
equation. We utilize the Markov Chain Monte Carlo (MCMC) ^ method to con- 
strain parameters of a model for multi-wavelength observations of TeV SNR RX 
J1713.7-3946. 



2. Stochastic Electron Acceleration in SNRs 



We consider the generic process of particle acc eler ation by compressible motions in 
SNRs as proposed by Bykov & Toptygin (1983)231. Schlickeiser (1984)1131 has shown 
that the evolution of the overall particle distribution function N can be studied by 
solving the diffusion-convection equation in the energy space: 
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where 7, -D 77 , T csc and Tuf e are the Lorentz factor of the particle, the diffusion 
coefficient in the energy space, the escape time from the acceleration region and 
the age of the remnant, respectively. Q is a source term and 7 corresponds to the 
radiative energy loss. Here 77 is related to the first-order Fermi acceleration and 
the terms related to Tiif e account for the effects of the remnant expansion with an 
expansion rate of 1/Tnf e . The evolution of N is then determined by -D 77 , 77, T csc , 7, 
and Q. 
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For the sake of simplicity, we will assume that Q and r\ are independent of time 
t, and for the study of electron acceleration, we consider the energy loss due to 
synchrotron and inverse Comptonization of the background radiation EHE1 F r 
stochastic interactions of electrons with the background turbulence, both _D 77 and 
T csc are determined by the turbulence properties. We consider a homogeneous and 
isotropic turbulence with a Kolmogorov spectrum, then the characteristic length 
of the magnetic field in a high (3 plasma is given by L — Ci(va/U) 3 R, where R, 
U , and va are the radius of the remnant, the shock speed, and the Alfven speed, 
respectively ®. C\ is a dimensionless constant. The spatial diffusion coefficient of 
relativistic particles is then given by 

X = cL/3 + a 7 m e c 3 /(3eB) (2) 

where we have assumed super-Bohm diffusion (with the dimensionless constant 
a > 1) for electrons with a gyro-radius r g = ^ra^c? /eB greater than L, and c, m e , 
e, and B are the speed of light, electron mass and charge, and the magnetic field, 
respectively. Then we have 

D 77 = j 2 U 2 /9X, (3) 
T CSC = R 2 /X. (4) 
To have a self-similar solution, T csc needs to scale with the age of the remnant T\a c . 

2/5 —3 /5 

For the Sedov solution of SNR evolution with R oc TJL and U oc T lifc ' , this implies 
that X oc RU, which, in combination with equation ([2]). gives va oc t/ 4 / 3 oc T,.^ 5 
and B oc T M 4/5 . 

For the shell Type TeV SNR RX J1713.7-3946, we adopt T ufe = 1600 
yr, R = 10(T lifo /1600yr) 2 / 5 pc, U = 3500(T lifc /1600yr)- 3 / 5 km/s, B = 
12.0(T lifo /1600yr)- 4 /VG, and assume that the electron acceleration starts at 70 
years after the supernova explosion L^- 1 and the electrons are injected at the proton 
rest mass energy with 7 = 1836. Kinetic plasma waves can acceler ate electrons to 
the proton rest mass energy efficiently at even lower energies HHHIHl Then by fitting 

1 /3 

the emission spectrum (Fig. [1]), we have r/ = —0.67, a = 1.3, C x va = 0.137J7, 
and the total energy of electrons with E to \ = 5.2 x 10 47 ergs. The corresponding 
reduced \ 2 °f the best fit is 1.83. In Figure HJ the radio datum is from Acero et 
al. (2009) 1201, other data are from Fan et al. (2010) EH. The right panel of Figure 
[T] shows the model predicted spectral evolution. Since the injection rate is constant 

—4/5 

and the magnetic field is proportional to T lilc , the radio luminosity has a very 
weak dependence on time while the GeV luminosity increases roughly linearly with 
time. The spectral hardening in the TeV band in the early phase of the evolution is 
closely related to the electron acceleration process. We notice that the self-similar 
solution may not be valid after 10000 years since the magnetic field falls below 3/iG, 
which is considered a typical value of the interstellar medium. Figure [2] shows the 
probability distribution function of model parameters given by the MCMC method. 
It is interesting to note that the value of a of the best fit is very close to 1 and T] is 
negative, implying a first-order energy loss process. 
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Fig. 1. The left panel shows the SED of the self-similar solution proposed in this paper at 1600 
years; The right panel shows the evolution of SED: 500 yr (dashed line), 1000 yr (dotted-dashed 
line), 1600 yr (thick line), 3200 yr (thick dotted-dashed line), 6400 yr (thick dotted line) and 12800 
yr (thick dashed line). 




Fig. 2. Probability distribution function of model parameters normalized at the peak value. 

3. Summary and Conclusions 

The F-P equation in general does not have analytical solutions for the diffusion coef- 
ficient in the SA theory. To obtain the electron distribution from the F-P equation, 
one may resort to some simplifying approximations involving only the time and 
energy variables. In the paper, we develop a numerical code to solve this reduced 
F-P equation, which is applicable for arbitrary forms of the diffusion coefficient. 

Recent observations of SNR RX J1713. 7-3946 made with the Fermi space tele- 
scope show that the spectrum of this source in the GeV band is very hard with 
a power-law photon index of T = 1.5 ± 0.1, which is difficult to accommodate in 
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a hadro nic scenario, but agrees well with the IC origin of 7-rays in the leptonic 
1991 

scenario 1 ^. On the other hand, by considering the potentially high inhomogeneity 
of the shocked interstellar medium (ISM), 1221 ar g ue that a hadronic model may 
still explain the Fermi observation. The nature of the TeV emission from SNR is 
an issue still open to investigation. In the work, we show that the SA of electrons 
by turbulent pla sma waves in the shock downstream might naturally explain these 
observations 1=^1 Future observations with the HXMT and NuSTAR can test the 
model. 
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